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Abstract 

We use a simple nonlinear scaler displacement model to calculate the distribution of effect created 
by a shear stress on a double stranded DNA (dsDNA) molecule and the value of shear force 
which is required to separate the two strands of a molecule at a given temperature. It is shown 
that for molecules of base pairs less than 21, the entire single strand moves in the direction of 
applied force whereas for molecules having base pairs more than 21, part of the strand moves in 
the opposite direction under the influence of force acting on the other strand. This result as well as 
the calculated values of as a function of length of dsDNA molecules are in very good agreement 
with the experimental values of Hatch et al. (Phys. Rev. E 78, 011920 (2008)). 

PACS numbers: 87.15.-v,64.70.qd,05.90.+m,82.37.Rs 
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I. INTRODUCTION 



A double stranded DNA (dsDNA) molecule consists of two polynucleotide strands con- 
nected loosely by hydrogen bonds through the base pairs and base-stacking between nearest 
neighbour pairs of base pairs and wound around each other to make a helix. The constraints 
of this helical structure require that the two base sequences on opposite strands must be 
complementary, with Adenine (A) always binding to Thymine (T) and Guanine (G) binding 
to Cytosine (C) [1]. The force that holds the complementary strands of DNA together is an 
important regulator of life's processes because the binding of regulatory proteins to DNA 
often involves the procedure of mechanical separation of its strands. The intermolecular 
forces of DNA have been studied extensively using a variety of techniques, which cover a 
broad range of forces from a few piconewton(pN) up to several hundred piconewtons (pNs). 

These techniques use either atomic force microscopy (AFM) or laser optical traps and 
magnetic tweezers {2 7|. The experiments fall into two general categories; those conducted 
on short DNA and those on long DNA. Experiments on long DNA focussed on the overall 
properties of the molecule js], Q] and resulted in the discovery of S-DNA js], loijlll: the 
occurance of this so called B-S transition has also been found in a short DNA chain of 
30 base pairs Those on short duplexes focussed on the reaction pathway of melting, 
and contributed to the understanding of the local unbinding of DNA 13l-ll7l|. In the single 
molecule experiments, the result may depend on choice of which variables are fixed and 
which can fluctuate. In some of the experiments (like AFM) extension is fixed and force is 
allowed to fluctuate (constant extension ensemble) while in others (like magnetic tweezers) 
the force is fixed and extension is allowed to fluctuate (constant force ensemble). 

When a force is applied to pull apart the two strands from one end of a dsDNA molecule 
in a direction perpendicular to the helical axis, bases are sequentially stretched as the duplex 
is unzipped. On the other hand, in shear unzipping in which the applied force pulls the two 
strands from opposite ends as shown in Fig.l, the stretching is spread out over many base 
pairs. The study of distribution of shearing force along the length of a DNA molecule may 
lead to valuable informations about the force distribution across the phosphate backbone of a 
single strand of DNA versus the force distribution across the paired bases of complementory 
DNA. Such informations are of relevance not only to understanding the biological processes 
but also in material science application such as determining the strength of DNA/gold 
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nanoparticle assemblies 



Following the work of Lee et al. 13| many workers 15 



16 



19 



- |26| measured the value of 



shear force which separates the two strands of a dsDNA molecule. It was found that the force 
increases linearly with the length. Using a simple ladder model of DNA and expressing the 
backbone bond energy as well as the interaction energy among complementary bases in the 
form of harmonic springs, de Gennes in 2001 predicted that the critical force for shear 
unzipping which for short chains shows linear dependence would saturate to a finite value 
in the limit where the number of base pairs approaches infinity and that the shear stress 
relaxes over a distance (number of base pairs) (~\/^' where k is the spring constant 
characteristic of stretching the backbone, and R is the spring constant characterstic of 
stretching the hydrogen bonds between base pairs) on either side of the chain. The critical 
force Fc for shear unzipping as a function of number of bound base pairs was given j27| as 



F, = 2Ux-Hanh{x-) + l]. 



(1) 



where fc is the rupture force of a single bond and (N+1) is the number of base pairs in the 
molecule. 

Chakrabarti and Nelson 18|] have generalized the simple harmonic model of deGennes 



271] by representing the interaction among complimentary bases by the Lennard- Jones (12- 



6) potential model and found that the strain is indeed localized over a narrow range of 
on either side of the chain and the chain unzips when the force exceeds a critical value of 
Fc ~ fcN for short chains and F^. ~ 2fcX~^ for long ones as predicted by deGennes ^^l^. It 
may, however, be noted that the calculations of deGennes 271] and also of Chakrabarti and 
Nelson Q correspond to .ero temperature. The Langevin dynamics srmulafon has recently 
been used to study shear unzipping of dsDNA at finite temperature [28| . The results found 



for Fc are in agreement with the results of deGennes [27| and Chakrabarti and Nelson 18]. 

Hatch et al. [26] have measured the value of Fc for several dsDNA molecules of length 
ranging from 12 to 50 base pairs at the room temperature and found that Fc is linear function 
of molecular length only up to 20 base pairs and approaches to an asymptotic value as the 
number of base pairs increases. In fact they found that Fc for a molecule of 32 base pairs 
has already reached within 5% of the asymptotic value which was found to be 61.4 pN. But 



when Hatch et al. 



261 ] tried to fit their data to Eq.(l) they found that this can be done only 



by assuming 7 base pairs of all chains to be in the open state, i.e. the effective length of a 
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molecule of (N+1) base pairs is (N-6) irrespective of the value of N. They attributed this to 
temperature effect as deGennes calculation did not include temperature. 

In this paper we calculate the value of F^. as a function of the number of base pairs at 



room temperature and compare our results with the experimental data of Hatch et al. [26 1 
and show that the need to adjust the length of a molecule is not as much due to temperature 
as due to use of the values of /c and x Eq.(l) and the nature of the curve of F^. vs N. 
We derive another form of Eq.(l) and show that with reasonable choice of values of /c and 
X one gets values of F^ which are in better agreement with experiment without adjusting 
the length of molecules than that found from Eq.(l). The model which we describe in Sec. 



II is similar to the nonlinear scaler displacement model of Chakarbarti and Nelson [18||. We 
calculate the force in both the constant force and constant extension ensembles. The paper 
is organized as follows. In Sec II we describe the model and calculational procedures. In Sec 

III we give results found for distribution of effects created by shear force along the length of 
a molecule and the value of Fc as a function of number of base pairs. In Sec IV we compare 
our results with those of deGennes 27| and with the experimental values 26|. The paper 



concludes with a brief comment given at the end of Sec IV. 
II. MODEL 

We consider a dsDNA molecule of length (N+1) base pairs which both 5'-ends (or both 
3'-ends) are pulled along the helical (molecular) axis by a force F as shown in Fig.l. The 
displacements of i^^ nucleotide from its equilibrium position are denoted by Ui for one (the 
lower one in Fig.l) single stranded DNA (ssDNA) chain and by Vi for the other (the upper 
one) chain. The effective Hamiltonian of the system can be written as 



N/2-l ^ N/2 

H= ^K[(ui+i - Ui)2 + (vi+1 - Vi)2] + v{\\i,-^,\) 

i=-N/2 i=-N/2 

-F.(uAr/2 - V_iv/2). (2) 

The first term of this equation represents the stretching energy of nucleotides, excecuting 
a simple harmonic motion with a spring constant k along each chains in dsDNA. In the 
second term, f(|ui — Vi|) represents the potential energy of interaction between bases in 
the i^^ pair and the last term of Eq.(2) represents the energy contribution due to the shear 



stress. In writing Eq.(2) we did not include the contribution arising due to the hehcity of 



DNA as this effect was found to be neghgible by Hatch et aL 26| and also by Lavery and 



Lebrun 



29|. 



For the potential f (|ui — Vi|) we use a simple model which has a hard-core repulsion and 
a long range attraction, 

7;(|ui-Vi|) = cx) for ^ < 0, 

for ^ > 0. (3) 



Here |zi| is the magnitude of increase in length of hydrogen bonds connecting bases in the 
z*'* pair from its equilibrium value, e is the depth of potential at the equilibrium separation 
and a is the diameter of dsDNA (see Fig. 2). The repulsion represents the steric hindrance 
which forbids the molecule from getting compressed along the bond linking the bases with 
respect to its equilibrium value. From Fig. 2 one finds , 

We take the value of a equal to 20A which is the diameter of the Canonical DNA (B-DNA) 
at room temperature. 

Since in shear unzipping, a molecule is stretched along its axis we consider the longitu- 
dinal displacements of nucleotides and neglect the transverse displacements; the transverse 



displacements have been found in ref 18| an order of magnitude smaller than those in the 



direction of shear. We define new variables, 

where, Ui and Vi now represent the longitudinal displacements of respective strands. In 
the notations used here Ui and Vi are positive when the i^^ nucleotide of upper and lower 
strands move to r.h.s. and negative when they move to l.h.s.. For the experimental situation 
shown in Fig. 1 it is clear that for i > 0, > and magnitute of Ui is greater than that 
of Vi and for i < 0, f j < and the magnitute of Vi is greater than that of Wj. This leads to 
following relation for variables Xi and Ui] 

x-i = -Xi, Xo = and y_i = yi. (6) 
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When we substitute these variables in Eq.(2) it decouples into two independent compo- 
nents; 



where 



H = H, + Hy, (7) 



N/2-1 ^ p 

o/^la^i+i - 7^(^Af/2 - X-N/2), (8) 

i=-N/2 ^ V 2 



and 



iV/2-1 ^ 7V/2 p 

Hy= J2 T^f^iVi+l - Vif + ^iV^) - -7^iyN/2 + y-N/2)- (9) 

i=-N/2'^ i=-N/2 

Here the potential f (|ui — Vi|) defined in Eq.(3) is expressed in terms of variable Using 
Eqs. (4) and (5), one can rewrite Eq.(3) as 



v(y.) = (10) 

Note that the expression of does not contain the on-site potential v{yi) and simply 
corresponds to a harmonic chain which is being pulled at the two ends by a force F/ a/2 
whereas the expression of Hy contains the on-site potential v{yi) as well as the force term. 

In view of the relations given by Eq.(6) the average value of displacement < x„ > of n*'* 
base pair can be calculated from the relation, 



<Xn>— N/2 , / N ' y^^) 



where /3 = (ksT)'^, being the Boltzmann constant and T is the temperature. Since for 
Hx given by Eq.(7) the integrals in Eq.(ll) are Gaussians, one can solve them analytically 
to give 

< Xn >= ^=-n. (12) 
v2k 



This result agrees with the one found by deGennes 27 1. 

The average value of displacement of n^^ base pair can be found from the relation, 



/ niI-jv/2 dViVn exp{-/3Hy) 



< yn >= . N/2 — : — :t7— • (13) 



The integrals appearing in this expression cannot be evahiated analytically because of the 
form of on-site potential v{yi). However, for the expression of Hy given by Eq.(9) the integral 
appearing in Eq.(13) reduces to multiplication of {N -\- 1) matrices. The discretization of 
the coordinate variable and introduction of a proper cut-off on the maximum values of y's 
determine the size of the matrices. We have taken -40A and 40A as the lower and upper 
limit of integration for each co-ordinate variable and discretized space using the Gaussian- 
Legendre method with the number of grid points equal to 900. By changing the limits of 
integration as well as the number of grid points we made it sure that the values of < y„ > 
are independent of the limit of integration and the number of grid points chosen to discretize 
coordinate variable. 

III. RESULTS 

When a shear force is applied on a dsDNA molecule, its two ssDNA strands get pulled in 
opposite directions as shown in Fig.l. The bonds in the backbone of DNA as well as bonds 
connecting bases in a pair are stretched. In a case of k being infinitely large the two strands 
will move like a rigid body pulling all base pairs in the sequence in parallel. The effect of the 
shear force will then be uniformly distributed across all of the base pairs. However, if k is 
finite, then both the backbone and the base pairs will stretch when a shear force is applied. 
The effect of the shear force may then be confined to limited lengths on both ends of the 
molecule. To see how the effect caused by shearing of a dsDNA molecule is distributed along 
the length of a molecule and how this depends on the energies associated with the stretching 
of backbone and base pairs we calculate the value of < y„ > from Eq.(13) for different values 
of n when the two end base pairs are stretched to a given length by the shear force. 

In Fig. 3 we plot our results for four dsDNA molecules of length 17, 25, 33 and 49 base 
pairs and for y-N/2 — yN/2 — 1.50A, 2.0A, 2.38A and 2.60A. The values shown in the figure 
correspond to e = 0.04eV, k = OAOeV/K^ and T = 300'^ K. When y-N/2 and yN/2 ai'c allowed 
to be free (i.e.F = 0), then < i/„ > is found to be zero for all n. From Fig.3 it is clear that 
for short chains the effect created by shear force is distributed along the entire length of a 
molecule affecting all base pairs, whereas, for relatively larger molecules the base pairs in 
the central part of molecules are only marginally affected. For example, for a molecule of 
length 17 base pairs the value of < yo > is 1.76A when y_8 — Vs — 2.38A, whereas, for 
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the similar situation (i.e. ?/_24 = 1/24 = 2.38A) the value of < ?/o > for a molecule of length 
49 base pairs is only 0.27A. The qualitative nature of these results are in agreement with 



the results found in refs [18| and [28|. The other point to be noted from the figure is that 
the qualitative nature of the distribution of the effect of shear stretching along length of a 
molecule is same for all values of y-.N/2 = yN/2 plotted in the figure. 

In Fig. 4 we plot the value of < ?/o > when y-N/2 = yN/2 = 2.38A as a function of length 
of dsDNA molecules. As discussed below, the value of y-N/2 = yN/2 = 2.38A is assumed to 
be the critical value of stretching in the sense that initiation of separation of two strands 
starts at this value of yN/2 and the shear force which creates this value of stretching of the 
end base pairs is equal to F^, the minimum ( or critical) force required to separate the two 
strands of a dsDNA molecule. From Figs. 3 and 4 it is clear that as one moves from either 
ends the differential force across the base pairs decreases and becomes very small at the 
centre for larger molecules but has not become zero at the centre even for a molecule of 
length 49 base pairs. The extension of the curve of Fig. 4 shows that it would be zero at the 
centre only for molecules of length larger than 80 base pairs. 

The average value of displacements < u„ > and < f„ > of ra*'' base pair can be found 
from the known values of < x„ > and <?/„>. From Eq.(5) one gets 

^ <Xn> + <yn> ,^ ^ <Xn> - <yn> /^,x 

< Un >= 1= ana < Vn >= 1= ■ (.14) 

V 2 V 2 

In Fig. 5 we plot the values of < > and < f „ > as a function of n for molecules of length 
17, 21 and 49 base pairs. The values plotted in this figure correspond to y-N/2 = yN^^ = 
2.38A and F=40.9, 47.4 and 60.6 pN , respectively for molecules of length 17, 21 and 49 
base pairs. These values of force, as is shown below, are the critical force of shearing of the 
respective molecules. From the figure we note that while for a molecule of 17 base pairs 
the entire ssDNA strand moves in the direction of applied force whereas in the case of a 
molecule of 49 base pairs nearly half of the strand moves in the opposite direction. While 
for a molecule of 17 base pairs, < us >= 0.66A, < Ms >= 2.70A and < V-8 >= 2.70A, 

< vs >= 0.66A, for a molecule of 49 base pairs < ^-24 >= — 2.86A, < ^24 >= 6.22A and 

< V-2A >= 6.22A, < U2A >= -2.89A. 

As long as the entire ssDNA strand moves in the direction of applied force the shear force. 



as explained in ref 18|], depends linearly on molecular length. The departure from the linear 
dependence of the shear force on molecular length is expected to take place when < U-n/2 > 
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and < f7v/2 > become zero i.e. the displacement of nucleotides on the opposite side of a 
ssDNA strand remains unaffected by the applied force. Indeed, we find that < m_io > and 
< viQ >=0 for a molecule of length 21 base pairs whcih is in very good agreement with the 
experimental value |26|. Since for molecules of base pairs larger than 21, part of a ssDNA 
moves in the opposite direction, a region develops inbetween on each strand which remains 
unaffected by the force. This region moves towards the centre on increasing the molecular 
length; the force gets saturated as soon as the region reaches the centre of the chain and 
stays there on further increasing the molecular length. 

The value of shear force needed to separate the two strands of a dsDNA molecule can be 
calculated in two different ways. In one, we follow a method proposed by deGennes [27| and 
which from hereon is referred to as a method of constant force ensemble. The other method 
is based on the constant extension ensemble. 

In the method of constant force ensemble one first defines a critical distance y for the 
rupture of a base pair (i.e. when ?/„ of n^^ base pair becomes larger than y the bases of 
the pair become free) and calculate the value of force which can stretch a base pair to the 
critical distance y. This force can be found from the on-site potential v{y). Thus 

If we take the value of y = 2.38A ,we find fc = A.lpN for the value of e and k, given above. 
This value of fc is close to the value used by Hatch et al (26| to fit their experimental data 
and the value estimated by Chakrabarti and Nelson [l8|. To rupture the end base pairs of 
a given molecule they must be stretched by the shear force to distance y. The balance of 
force at one ends of one of the ssDNA gives 



Fc = k{< Un/2 > - < Un/2-1 >) + fc- (16) 

Using the relations of Eq.(5) and Eq.(12) we get 

Fc = V2k{< yN/2 > - < yN/2-i >) + (17) 

Taking the value of < y-N/2> =< yN/2 >= y =2.38 A we calculate the values of < 
y N/2-1 > from Eq.(13) for several molecules of length 10-60 base pairs. The value of force 
Fc found from Eq.(17) is shown by dotted line in Fig.6. 
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In the constant extension ensemble one first calculate the work done in stretching the 
two end base pairs of a given molecule to distance y. This work is found from the relation 

W{y) = ^[InZ^+M - /nZjv+i], (18) 

where 

ZN+i{y) = n dy^6{y^N/2-y)S{yN/2-y)exp[-l3Hy], (19) 

i=-N/2 

is the constrained partition function and 



N/2 
[[ 

-N/2 

is the partition function. 6 is the Dirac function. 



Zn+i = j n dyi^M-I^Hy]. (20) 



We used the matrix multipliction method to evaluate W{y) from the above equations for 
e = 0.04el^, k = O.leV/A'^ and T = 300^ K. The derivative of W{y) with respect to y gives 
the average force F{y) that is needed to keep the extension of end base pairs of a given 
molecule equal to y. Thus 

To get the value of critical force we have to chose a value of y which corresponds to rupturing 
of base pairs. The values shown in Fig. 6 by dashed line are found when y was taken equal 
to 2. OA; this value is slightly lower than the one taken for the constant force ensemble. The 
difference in the value of critical stretching in the two ensembles may be due to difference 
in the path of unzipping. 

We note that the values of Fc found by methods of the constant force and the constant 
extension ensembles are close but not identical. The difference between the two as expected 
30l | is large for molecules of smaller lengths but they becomes close as molecular length 
increases. Both methods give the same asymptotic value, equal to 61.2pN, which is in very 
good agreement with the experimental value 61.4pN [26|. The experimental values shown 
in the figure are of Hatch et al 26|. In view of large spread in experimental data we find 
good agreement between experimental values of and the theoretical values found using 
the constant force ensemble as well as the constant extension ensemble. 
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IV DISCUSSIONS 



(a) Comparison with the experiment 

Since the experimental values 



26| of Fc given in Fig. 6 are found using the fixed force 
ensemble, we concentrate our discussion with the theoretical values found using the same 
ensemble and shown in Fig. 6 by dotted line. The agreement between theory and experi- 
ment is excellent except for a molecule of length 12 base pairs. We know that at a given 
temperature shorter dsDNA molecules are less stable compared to longer molecules. It is 
quite possible that at room temperature due to surface effects few of base pairs at the two 
ends of a molecule of length 12 base pairs are nearly in open state. The curve shown in 
Fig.6 by solid line is found when the effective length of a molecule of length 12 base pairs 
is taken to be equal to 10 base pairs and that of a molecule of length 16 base pairs equal 
to 15 base pairs. The values of force Fc shown by solid line is in excellent agreement with 
;he experimental values for the entire range of molecular length investigated by Hatch et al. 



26| . The saturation value of F^. found to be 61.2 pN is also in very good agreement with 



the experimental value of 61.4 pN [26]. 



(b) Comparison with the results of deGennes 27|] 



In order to compare our results with those of deGennes [27] we first estimate the value 
of X defined as x = \f^i where R is the spring constant of a simple harmonic potential 
between the bases of a pair. Expanding v{y) in ascending powers of y one gets 

<m) = -^AC^)yl (22) 



2 Va2. 
„ 12e ,2 24e 

R = — and X = (23 

Substituting the values of e, n and a given above we find = 6.44 which corresponds to 
k/R = 83.3. This value of k/R is in good agreement with the predicted value of 77 based 
on calculation of the spring constants for base pairs and backbones 13, 

The expression for the displacement < ?/„ > found by deGennes |27| can be written as 

<yn>=l<yo>{e^'' + e-n- (24) 
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Taking the value of x^^=6.44 and the values of < ?/o > determined above (shown in Fig 4) 
we calculate the value of < ?/„ > from Eq.(24) for molecules of length 23 and 49 base pairs 
and compare them in Fig 7 with the values of < y„ > found from Eq.(13). The two values 
do not agree; the difference increases with the length of molecules. A good agreement is, 
however, found when x~^=9.4 is taken. This value of is about \/2 larger than the value 
found above. We calculate Fc from the deGennes equation given by Eq.(l) taking x~^=6.44 
and 9.4 and /c=4.1 pN. The results are plotted in Fig 8 in which we also plot experimental 
values and the values found from our approach and shown in Fig 6 by dotted line. The 
values found from Eq.(l) are very different from both experimental values and values shown 
by dotted line. In ref 28| the value of Fc/ fc as a function of N were found to be in good 
agreement with the values found from Eq.(l) for = 10. 

Since Eq.(l) involves the rupture force fc for a base pair and x, the value of Fc depends on 



the values of these quantities. Most often quoted value of fc is in the range of 4-5 pN |18 



and that of x ~ 0.1 [l8|, l26|-|28|. Hatch et al. [26] took /c=3.9 pN and x"^=6.8 to fit their 
data to Eq.(l) and found that they have to adjust the length of molecules which amounted 
to shifting the curve of Fc to the right direction by 7 unit of base pairs in order to get good 
agreement with the experimental values. We now examine whether such an adjustment of 
length is essential or there is some other reason for not getting good agreement. 

In the appendix we derive Eq.(24) from Eq.(9) and show that x which appears in Eq.(24) 
should be defined as X=\[^- Substituting the values of k and R we find x^^=9.2 which is 
close to 9.4 found by using the values < y„ > found from Eq.(13). With this definition of x 
we find the following expression for Fc (see Eq.(A8)). 

Fc = V2fc[x'Hanh{x^) + 1]. (25) 

The appearance of a/2/c instead of 2/^ (in Eq.(l)) is due to difference in the definition 
of X- The definition x = \[^ seems more appropriate in the sense that it is ratio of the 
two harmonic spring constants than the definition x = \/v H which one spring 

constant is multiplied by two. 

The value of Fc found from Eq.(25) when x~^=9.4 and /c=4.1 pN are compared in Fig 9 
with the values using the method of Sec. Ill and shown in Fig 6 by dotted line and values 
which are found from Eq.(l) when x^^=6.8 and /c=3.9 pN j26|. From the figure we first 
note that Eq.(25) gives values of Fc as a function of molecular length which are close to the 
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values found in Sec. Ill (dotted line); the difference between the two is due to the combined 
effect of nonlinearity and the temperature. Thus the effect of temperature is not as large 
as has been suggested in ref [26|. The necessity for adjusting the length of molecules is due 
to the values of x ^-nd fc used in Eq.(l) by Hatch et al. |26|. If one takes x~^=9-4 and 
/c=4.1/-\/2 pN and calculate Fc using Eq.(l) one gets values shown by full line in Fig 9 
which are found from Eq.(25) with x~^=9.4 and /c=4.1 pN. 

In conclusion, we developed a method to calculate the distribution of shear force along 
the length of a dsDNA molecule at a finite temeprature. The value of shear force which is 
required to separate the two strands of a molecule has been calculated in both the constant 
force and the constant extension ensembles. The values of Fc found by these methods 
differ for molecules of shorter length but approach to each other as length increases. Both 
methods gave the saturation value equal to 61.2 pN which is in very good agreement with 
the experimental value 61.4pN. The value of Fc is found to increase linearly with length up 
to 21 base pairs in agreement with experimental results. The plots of < u„ > and < f „ > 
given in Fig. 5 show that as long as the applied force pulls the entire strand in its direction of 
application, Fc depends linearly on length of molecules. The departure from linear behaviour 
takes place when part of a strand moves in opposite direction under the influence of force 
pulling the other strand. The saturation value is achieved when half of a strand moves in the 
direction of force applied on it and the other half in the opposite direction. It is shown that 
the value of Fc quickly attains its saturation value on increasing length of DNA molecule; 
for a molecule of 32 base pairs Fc= 57.0pN as compared to saturation value 61.2pN. 

The agreement between theorectical and experimental values of Fc shown in Fig. 6 indi- 
cates that the model proposed in this paper is capable of describing the responses of dsDNA 
molecules to shear stress. The model, however, neglects the effect of helicity of DNA and 
has assumed the molecule to be homogeneous. As far as the effect of helicity is concerned, 



it has already been shown to be negligible by Hatch et al. [26j and by Lavery and Lebrun 
2^. The effect of heterogenity arising due to random distribution of A-T and C-G base 
pairs in a sequence is also expected to be small as long as the shear unzipping from the two 



ends is symmetric, 
gated by Hatch et al 



'his is corroborated by the fact that although the molecules investi- 



261] are heterogeneous with half G-C and half A-T base pairs, yet their 



response to shear stress is described very well by a model which assumes molecules to be 
homogeneous. The values of parameters e and k will, however, depend on the percentage of 
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A-T and C-G base pairs in a given sequence and on the temperature. But, in the absence 
of symmetry in the unzipping from the two ends a quahtative new features, as argued in ref 



181], may arise. 
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Appendix 



In this appendix we derive an expression for tjn from Eq.(9). Writing f (?/„) as (see Eqs.(22) 
and (23)) 



v{yn) = -e + -Ryl 



(Al) 



where R=^ and differentiating with for — ^ < n < ^ we get the following equilibrium 



conditions 



27| from Eq.(9) 



dH 



dyn 

For n = N/2, one gets 

dH ^ ^ F 

From these equations we find 



{A3) 



where x=\i^^ and 



F 1 

^° " V2Rcosh{lxN)[x-'tanh{lxN) + 1] ^^^^ 

Let the force on the last hydrogen bonds {n=N/2) is RyN/2 when we reach the threshold fc 
( fc being the rupture force for a base pair). Thus 

/, = Ryocoshi^xN) (AQ) 

From Eq.(A3) we get 

Fc = V2KxyN/2tanh{lxN) + v^/c {A7) 



V2ax-Hanh{^xN) + l] {Al 
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This equation differs from Eq.(l) in the definition of % and the mutiplying factor which is 
now \f2 instead of 2. 
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FIG. 1: Schematic of a dsDNA molecule of N+1 base pairs under a shear stress. Un and Vn are 
displacements of n^^ nucleotide in lower and upper strands, respectively. 
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FIG. 3: Displacement < y„ > as function of n for four molecules of length 17, 25, 33 and 49 base 
pairs and for y^N/i = yN/2 =1-5A (full line), 2.0A (dashed line),2.38A (dotted line) and 2.6A 
(dot-dashed line). 




FIG. 2: Stretching of a base pair under the influence of shear stress. The open circles in (a) indicate 
the equlibrium position and filled circles stretched positions. 
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FIG. 4: Displacement < yo > of the central base pair as a function of number of base pairs of 
dsDNA molecules. 




n 



FIG. 5: The displacements < > and < u„ > as a function of n of molecules of 17, 21 and 49 
base pairs for the shear force Fc= 40.9, 47.4 and 60.6 pN, respectively. 
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FIG. 6: The value of critical force in pN as a function of number of base pairs of a dsDNA 
molecule. The curve shown by dotted line corresponds to constant force ensemble and dashed line 
to constant extension ensemble. The curve shown by full line corresponds to values found by taking 
the effective length 10 and 15 base pairs for molecules of length 12 and 16 base pairs respectively. 



The experimental values of Hatch et al 



26l | are shown by diamond when 5 '-ends are pulled and by 



star when 3'-ends are pulled by shear force. 
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FIG. 7: Comparison of the values of < y„ > found from Eq.(12) (full line) and the deGennes relation 
Eq.(24) (dashed line) when x~^=6.44. For x~^=9.34, the values found from Eq.(24) overlap with 
the curve drawn by full line. 
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FIG. 8: Comparison of the values of shear force as a function of the number of base pairs in 
a molecule. The curves shown by dashed and dot-dashed lines represent the values found from 
Eq.(25) when x~^=9.34 and 6.44, respectively. The dotted line represents corresponding curve of 
Fig.6. 
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FIG. 9: Comparison of the value of shear force as a function of the number of base pairs in a 
molecule. The curve shown by dot-dashed line is found from Eq.(l) with x~^=6.8 and /c=3.9pN 
26[, the curve shown by full line is found from Eq.(25) with x^^=9.4 and /c=4.1pN and the dotted 
line represents corresponding curve of Fig. 6. 
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